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Abstract 

In an effective operator approach, the full set of leading contributions to anomalous top couplings 
comprises various new trilinear as well as higher interaction vertices, some of which are related to 
one another by gauge symmetry or equations of motion. In order to study trilinear top couplings 
to SM gauge bosons such as tt'j, ttZ, tbW and ttg, the operator set can be restricted accordingly. 
However, the complete basis cannot be consistently mapped onto an on-shell parametrisation of the 
trilinear vertices alone. Four-fermion contact terms qqtt and udtb must be included if consistency 
is to be retained. In this paper, we point out how these interactions contribute to the single top 
search channels for anomalous trilinear tbW couplings at the LHC and Tevatron, thus affecting the 
corresponding bounds. All results are based on full leading order partonic matrix elements, thus 
automatically accounting for off-shell and interference effects as well as irreducible backgrounds. A 
discussion of the quantitative effects of going from on-shell to off-shell tops and applying acceptance 
cuts on the matrix elements is also provided. 
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I. INTRODUCTION 



The Standard Model of Particle Physics (SM) has stood its ground during the past 
decades with great success, consistently explaining and predicting a great variety of high 
energy experiments with unchallenged precision. One of the major cornerstones was the 
discovery of the top quark at the Tevatron in 1995 [1, 2], confirming the postulated three- 
family doublet structure of the SM. While the Tevatron experiments have continued to collect 
data and improve their measurements of top properties, most importantly its mass [3- 
6], most attention is now directed to the LHC up and running at y/s = 8TeV, and the 
results of its multi-purpose experiments ATLAS and CMS improving on top statistics by the 
day. By now, top pair production has been measured in different channels with remarkable 
accuracy [7-12], and also single top production has already been established for bg — > tW 
associated production despite its small cross section [13, 14] resp. definitely observed in the 
dominant t channel bq — > tq' [15-17]. The ever growing abundance of top events at the LHC 
is beginning to allow the determination of more involved observables such as asymmetries, 
invariant pair mass distribution and top couplings to the other SM particles with high 
precision (cf. e. g. [18, 19] for an overview). 

On the theoretical side, the top quark takes an outstanding place among the spectrum of 
SM particles as a possible window to new non-SM physics because of its uniquely large mass 
of the order of the electroweak symmetry breaking (EWSB) scale m t ~ v ~ (9(100 GeV), 
with its role within the dynamics of EWSB still unresolved. Corresponding new physics 
effects in the top sector may manifest themselves in the deviation of top properties from 
their SM values, where the main attention in this paper is directed towards the trilinear 
couplings to SM gauge bosons, especially the charged-current (CC) interaction tbW. There- 
fore, even before any experimental analysis, a theoretically robust parametrisation of these 
anomalous couplings has to be found, at the same time reducing the parameter space to an 
experimentally manageable minimum while staying fully general within the basis of effective 
operators generating these couplings at Lagrangian level. Indeed, starting from the complete 
set of effective dimension six operators as written down by Buchmiiller and Wyler already 
in 1985 [20], substantial effort has been put into this task in the past decades by various 
authors [21-30]. The crucial ingredient of most of these analyses is to employ the theo- 
rem [31-35] that the field equations of motion (EOM) can safely be used at the Lagrangian 
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level at a fixed order in the effective operator expansion in order to rewrite operators and 
identify redundant structures. 

As a result of this procedure, it is often argued that these redundancies allow for a 
reduction of independent couplings to be incorporated in a phenomenological analysis of 
anomalous top couplings. However, as has been pointed out e. g. in [27, 28] and will also 
be reviewed again in more detail later on, the application of the EOM necessarily generates 
four-fermion contact interactions, which are nevertheless often dropped from the analyses for 
the sake of simplicity. We note that the latter procedure does not correspond to a rewriting 
but rather to a redefinition of the originally chosen operator basis, thus also departing 
from the full generality of the original basis with respect to the richness of structures in 
the trilinear couplings. Still, the operator equalities derived by various authors [27, 28, 36] 
and systematically presented in [28] are very useful to simplify an implementation of the 
most general set of trilinear top couplings into a Monte Carlo (MC) generator in a gauge 
invariant way, so in our approach, rather than dropping part of the physics, we make use 
of these equalities to implement all the trilinear top couplings to SM gauge bosons in the 
language of on-shell couplings including the required quartic contact terms into the parton- 
level MC event generator Whizard [37], also addressing the interplay of anomalous top and 
bottom couplings — the latter already heavily bounded by LEP data — and the repercussions 
on the top couplings. Finally, we present phenomenological consequences obtained with our 
implementation for the parameter space of the anomalous couplings in the CC sector. 

This article is organized as follows: in section II we review the procedure described 
above of defining a complete operator basis to generate anomalous top couplings to SM 
gauge bosons and applying the EOM to rewrite some of these operators, thus arriving 
at the most suitable form for a MC implementation. In section III we discuss the LHC 
phenomenology with a focus on single top production, including a comparison of on-shell 
and off-shell approaches to retrieve the cross sections at detector level as functions of the 
anomalous tbW couplings as well as a presentation of the physical effects and consequences 
of the newly added coupling structures. A discussion and summary of the main statements 
and results can be found in section IV. 
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II. THEORETICAL SETUP 



In order to be self-contained in this article, we start this section by reviewing in some 
detail the main steps of the procedure presented in [28] to simplify the most general set of 
operators generating the trilinear anomalous top couplings to SM gauge bosons ttj, ttZ, 
tbW and ttg. However, although we emphasize that all of these couplings are implemented 
in Whizard in a gauge-invariant way, including all quartic terms generated by the operator 
rewriting, we will restrict ourselves here to the discussion of the CC sector, i. e. only those 
operators generating anomalous contributions to the tbW interaction, for two simple reasons: 

1. the complexity of the parameter space in the CC sector is increased in a minimal way 
compared to previous studies [36, 38-40], because it turns out that only one additional 
operator (and hence coupling) has to be considered; 

2. at hadron colliders, the experimental access to the new effects is most straightforward, 
combining studies of CC single top production and of top decay products 1 . 

In the following two subsections we develop the basic ingredients of the effective operator 
analysis, recapitulate the operator rewriting procedure and finally present our extended 
parameter space for the anomalous tbW couplings. 

A. Effective Operator Approach and Operator Basis 

There are basically two ways to tackle new physics beyond the SM in a systematic and 
consistent manner: Either the model building (top-down) approach, i. e. starting from a 
postulated Lagrangian — which incorporates a sensible UV completion — and deriving from it 
physical effects to which present or planned experiments might be sensitive, or the effective 
(bottom-up) approach, i. e. starting from the established SM symmetries and a priori con- 
sidering all possible new physics effects compatible with these symmetries at the Lagrangian 

1 The anomalous NC sector, while of course related to the CC sector by gauge symmetry (cf. end of 
Sec. II B), is much harder to access experimentally, because one would have to identify the final state ttZ, 
which is an even more complex analysis than the already challenging tt-f study due to the further reduced 
cross section and the necessity to reconstruct the decaying Z. In the QCD sector the deviations from the 
usual parametrisations as studied e. g. in [41-43] appear in the qq — > tt amplitudes only and are therefore 
suppressed by the pdf 's with respect to the dominant gluon fusion channel. 
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level, postponing the question which larger theory might generate the relevant parameters 
at a higher energy scale A. 

Since we follow the second approach, it shall be clarified a little further. The idea is 
to confront new physics completely unbiased, that is without any assumptions about the 
dynamical degrees of freedom generating it, and to study the effects that are manifest at 
a testable energy scale (considerably smaller than the resonant scale A) where the degrees 
of freedom are the well known SM particles. This corresponds to integrating out the heavy 
modes, thus generating effective operators of mass dimension d > 4 which are nor- 
malized by appropriate powers of A. In the model-independent approach, the effective 
Lagrangian can be written as an expansion in 1/A [20, 44, 45]: 

r (d) 

d>4,i 

with dimensionless operator coefficients C\ , comprising all possible effective operators built 
from SM fields and derivatives only, and compatible with all local and global SM symmetries. 
A complete set of these operators for d = 5, 6 can be found in [20]. 

The only possible d — 5 operator in this setup is a neutrino mass term [20], so the leading 
contributions to anomalous trilinear top couplings must be d = 6. The complete operator 
list at this order can be found in [28-30] , of which we now quote the ones relevant to trilinear 
tbW interactions (also adopting the nomenclature of [28]): 

Off = i (0 V/^) {qjrf 7*014) , (2a) 

0% = i^D,cf){u R ^d Rj ) , (2b) 

Oii w ={q L ^T I u R] )4>W I u , (2c) 

OH w =(q Li a^7 I d Rj )cj ) Wi u , (2d) 

0% v =(q La '"7 I D^ qLj )W I u , (2e) 

with generation indices i,j = 1, 2, 3 and the non-Abelian SU(2)^ field strength components 

K» = W - d»K - g e IJK WtW v K (3) 

to be contracted with the Pauli matrices t 1 (I = 1,2,3). The qL(R)i are left (right )-handed 
quark spinors in the electroweak isodoublet (isosinglet) representation, and <fi is the isodou- 
blet complex SM scalar field acquiring a vev (0) = ^(0,"u) T , and = ir 2 (fi*. Of the 
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other electroweak operators listed in [28], and OV u as well as all those containing the 

hypercharge field strength only contribute to NC interactions, whereas the operators 
Op tt , , O^d an d % bd a PP ear to contribute to the tbW vertex. However, the differences 



W Du(d) 



^Du(d) 



are entirely redundant as is shown in [28] , and the sums O l ^ u ^ + 0^ u(d] are 



Du{d) 



proportional to the gauge boson momentum q^ = (pi — PjY so that amplitudes containing 
these vertices vanish either for physical on-shell W or for on-shell light fermions coupling to 
the W, which is always the case at parton level for all processes to be considered for single 
top and top decay studies discussed here. So Eq. (2) represents the most general d = 6 
operator basis generating anomalous tbW couplings, in which we shall therefore be complete 
in our phenomenological studies. 

We could now straightforwardly start off to find and implement all interactions generated 
by the operator basis which could appear in the relevant amplitudes. However, this can 
become a rather involved business particularly for O l ^ w which contains, apart from the 
trilinear coupling, also some relevant quartic terms such as e. g. tbWg with a complicated 
Dirac and momentum structure. Therefore, in order to facilitate the implementation work 
we rather follow the operator rewriting procedure of [28], illustrating the main steps here for 
O l q W : starting from its decomposition into hermitian and anti-hermitian parts, the hermitian 
part becomes 



(4) 



(dropping the total derivative), where the EOM of the W field 



(D v W u tf = g \ha 



(5) 



2 2 

can be applied to replace the derivative. On the other hand, with some algebra [20, 28] the 
anti-hermitian part can be brought in the form 

1 



(Olw) 
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(fca^ fcj f;-h.c.) 



(6) 



up to total derivatives, where the EOM of the quark field 



i&q L i = Y%u Rj t + Y*dn j <i> (7) 
(introducing Yukawa matrices Y u ^ d ) can be inserted. Joining it all together, one arrives at 
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the operator equality 



W qW 



+ | [gLil^r 1 q L j) (hkJi/r^Lk) + {iLkl '^r 1 ' q L k) 



+ 
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4 L 
1 

4 L 



O 



<f>q 



+ wry 



•yu r\ik _i_ yd r\ik 
1 jk w uW 1 jk w dW 



(8a) 
(8b) 
(8c) 



Obviously, the terms in (8b) and (8c) are redundant and can be absorbed into the opera- 
tors (2a)-(2d), whereas (8a) generates four-fermion contact interactions. 

Although it is clear that without further restrictions there is enough freedom within 
the operator basis to independently vary all the couplings emerging from (2a)-(2d) and the 
associated contact terms coming from the rewriting (8), the rewriting procedure corresponds 
to a shift of the original operator coefficients. Setting % = j = 3 and dropping all generation 
superscripts from now on, these shifts are: 



5ReC$ 



81m C uW 
81m C dW 



; ReC, 
m b 



V2' 



v 



qW , 

ImCgw , 
Im C qW ~ 



(9) 



assuming an approximate decoupling of the third generation in the Yukawa matrices. With 
this setup, we can now go on to physical states of the gauge and matter fields and write 
down the interaction terms generated by our operator basis. 



B. Parametrisation of Anomalous Charged-current Couplings 

Inserting the scalar vev and physical states of the gauge fields into the operators (2) 
and forming hermitian combinations C x O x + C*O x , one finds various trilinear interaction 
terms tbW, ttZ, ttA but also bbZ and bbA, as well as associated quartic interactions ttWW, 
bbWW, tbWZ and tbWA which are all necessary to maintain gauge invariance in the resulting 
amplitudes, and have therefore been included in our implementation. The resulting effective 
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tbW interaction Lagrangian can be written as 

Caw = ~ ^=b-f(V L P L + V R P R )tW~ + h.c. (10a) 

- ^ (g L P L + g R P R ) t W~ + h.C. (10b) 

~ ^^^(VE^tw- + h. c . , (ioc) 



where all couplings except for Vl = V t b — 1 vanish in the SM at tree level, and get the 
following anomalous contributions from operator coefficients 2 : 

SVl = (Cf> + 9 -Re C qW ) ^ , 6g L = V2C* dW ^ 2 , 

SV R = ^C* H — , 5g B = V2C uW — , 

5Vf= 9 -ReC qW ^ 2 . (11) 

The interaction terms (10a) and (10b) represent the on-shell parametrisation widely used 
in various phenomenological studies (normalization convention taken from [28]), which is 
retrieved from the operators (2a)-(2d). The interaction (10c) emerges from the hermitian 
part of O q Wi 

O q w + OqVK = (<?L37 m t 7 <?L3) (^W*) + higher contact interactions, (12) 

cf. Eq. (4), which — unlike the anti-hermitian part — cannot be completely recast into a 
combination of the other four operators. However, the partial redundance of O q w has been 
made explicit in the parametrisation (10) by defining its on-shell part into Vl so that any 
contribution ~ V£ s vanishes when the W goes on the mass shell. Hence it is no surprise 
that in 5Vl of Eq. (11) we find again the shift of the coefficient Ci^ already stated in 
Eq. (9) after the operator rewriting. Furthermore, by comparison to Eq. (8), one finds that 
all contributions ~ V£ s must be in one-to-one correspondence to the four-fermion contact 
interactions given in (8a), which is also highlighted by the fact that in physical amplitudes 
the kinematic structure of the W propagator is exactly cancelled by the g-dependent vertex. 

We have now isolated the non-redundant contribution of O q w to the tbW interaction 
Lagrangian, and also identified the most convenient way to implement it in a gauge-invariant 

2 Note that in Eq. (37) of [28] the operator coefficient C|| appearing in SVr should also be complex- 
conjugated. 
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way, namely by adding the quartic fermion vertices 



AC =i^(WP L t) [(u kJfl P L d k ) + (u k ^P L e k )] + h.c. (13) 
with g x = g Re C qW 

(cf. [30, 46, 47]), giving a relation of coefficients 

Of course, one might ask at this point if such a coupling structure should be counted among 
the anomalous tbW sector, but then again it must be noted that as a consequence of the 
common operator basis, the trilinear couplings are related to V£ s through the underlying 
operator coefficients. Specifically, Eq. (11) illustrates that a limit on 5Vl cannot be un- 

(3) 

ambiguously mapped onto a limit on the operator coefficient C\ without also bounding 
5V£ S (or the anomalous NC sector, see below). Moreover, the operator basis (2) and the 
corresponding set of couplings (10) parametrise all anomalous diagram insertions which can 
interfere with the SM diagram in a minimal way, making this approach consistent at the 
amplitude level. Finally, as pointed out in Sec. Ill C, the inclusion of the additional coupling 
also affects the interpretation of current and upcoming experimental results at the LHC. 

Before moving on to the phenomenological implications, let us discuss briefly the issue of 
anomalous bottom couplings within the effective theory approach: Since the original effective 
operators by construction respect the full electroweak gauge symmetry SU(2)^ x U(l)y, it 
is no surprise to find certain relations within the set of anomalous electroweak couplings of 
the heavy doublet (t, b) after spontaneous symmetry breaking. For example, an anomalous 
CC contribution 5Vl is directly related to the anomalous left-handed NC vector couplings 
ttZ and bbZ, the latter one stringently constrained by LEP data, so turning on 5Vl while 
respecting all existing bounds necessarily implies a non-vanishing anomalous contribution to 
the left-handed ttZ vector coupling [28, 48], or a fine-tuned relation with 8V£ , cf. Eq. (11). 
Similarly, Sg^ is directly related to the anomalous ttZ/ttj tensor couplings, just like Sg^ 
is to the bbZ /bbj ones (cf. e. g. [28, 49] for details). In short, it is impossible to vary 
the anomalous CC couplings in a consistent way within the effective operator approach 
without either getting anomalous NC couplings or including additional operators to fine- 
tune these effects away. Although these relations basically have no effect on a purely CC 
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FIG. 1: Diagrams contributing to LO on-shell single top production (anomalous tbW vertex 
marked by a dot): s channel tb production (left), t channel tj + tbj production (center) and 
associated tW production (right). 

single top study, one should bear them in mind when addressing anomalous CC couplings 
(the Whizard implementation contains the option to automatically enforce these relations). 

III. LHC PHENOMENOLOGY 

Apart from indirect searches using low energy observables, e. g. in the meson sector 
(cf. [50] for instance) or SM precision observables [51], there are basically two different classes 
of direct observables for top quark properties at the current collider experiments, namely 
those related to top production or top decays. While it is clear that only a combination 
of all available observables will deliver the best bounds on anomalous contributions, it is 
crucial to understand each analysis separately before the combination step. Therefore, we 
will focus here on the discussion of single top production cross sections, citing and using 
results from top decay studies to derive estimates for the most stringent bounds on the full 
anomalous parameter space at the end of the article. 

Single tops are produced at the LHC (and Tevatron) in three different channels, namely 
s channel tb production, t channel tj production (where j denotes a light hadronic jet), 
and associated tW production, cf. Fig. 1. While experimentalists are struggling to identify 
and discriminate these channels at the detector level with suitable selection criteria, the 
theoretical question is how the corresponding measured cross sections af ct (for final states 
% = tb, tj, tW) are represented as functions on the anomalous parameter space, i. e. how the 
measurement can be converted into bounds on the parameters. In this respect, a first step 
may be to separate the detector response from the hard production cross section: 




(15) 
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summing over partonic production processes j as functions of the parameter point g. de- 
notes the detector efficiency matrix mapping the process j onto the final state selection i, 
which can be retrieved with a detector simulation. Once the functions cxf art (g) are known, 
experimentally measured confidence intervals for af et can be mapped onto confidence inter- 
vals for g by formal inversion of Eq. (15). 

However, the remaining question to be addressed in this approach is: Where did we put 
the detector acceptance into e or into crf art ? The significance of this question is obvious, 
since anomalous couplings might very well affect the differential distributions, thus making $ 
a function of g. Therefore, the answer to that question influences the strategy as well as the 
efforts necessary to compute the function cr? art (<7), and potentially also the bounds derived 
from it, as we will show in the following section. 

A. Technical Setup 

On-shell Approach 

The simplest approach is to neglect the ^-dependence of the acceptance entirely and 
pull it into £, implying that the crf ar (g) in Eq. (15) represent the set of total partonic cross 
sections integrated over the full phase space. Further neglecting finite width and interference 
effects with irreducible backgrounds enables one to decompose of art (g) as 

rt=rt>n»». as) 

where of rod ((f) denotes the full on-shell single top production cross sections, and the product 
of branching ratios accounts for the decays of the heavy particles, namely t and one or two 
W's, depending on the production channel. Since all the on-shell production diagrams can 
contain only one anomalous tbW vertex insertion, it is argued in [36] that af rod (g) may be 
written as a polynomial up to second order in g: 

k,i 

where the of M are the total SM cross sections, and the K l kl denote the integrated products of 
diagrams with one insertion of and gi each, normalized to the SM point in each production 
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channel i. Plugging (16) and (17) into (15), one arrives at the ansatz employed in [36]: 



j,k,l 



■ 17 j M ■ «fci 9k 9i 



= Y, £ «-°r-4i9k9i, (18) 

j,k,l 

to be referred to as on-shell approach from here on. The advantage of the formula is obvious: 
once e and the k's are known, the conversion of measured results into bounds on g becomes 
very simple and efficient. However, the validity of the assumptions leading to this result 
shall be addressed now. 



Off-shell Approach 

Eq. (18) tells us that one should be able to vary the coupling point g within the ranges 
relevant for the study, with only minor effects on the detector response e in the phase space 
window which corresponds to a given final state selection. However, Eq. (18) implies even 
more, namely that retrieving the matrix element response as a function of g and applying 
acceptance cuts on the phase space should approximately commute, or equivalently, Eq. (18) 
should give the same results as e. g. 

°t ct {9) = E4 • [ $Part x (9) , (19) 

j 

where the basic detector acceptance cuts such as pt and rj cuts on the partons and leptons, 
represented by $ part , are applied directly to the phase space integration and hence formally 
included in the (/-dependent part of the formula, while the matrix e', assumed to be constant 
in g, denotes the efficiency of mapping the partonic final states at the acceptance level $ part 
onto the final state selections at detector level. 

To be more explicit, the idea is to accommodate as much of the acceptance cuts as 
possible within the (/-dependent part without becoming exclusive to any of the different final 
state selections, which are still contained in the ^"-independent e'. This obviously implies 
that the phase space window covered by e' must be fully contained within the acceptance 
window <3> part , leading to the notion that partonic acceptance and final state selection cuts 
should be adapted to each other as closely as possible. 
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Assuming leptonic t decay, we therefore apply the following acceptance cuts on the par- 
tonic phase space integration: 

$ part : p T (£, v) > 25 GeV and \rj (£) | < 3 , (20a) 

p T (j, b) > 30 GeV and \r}(j,b)\<5, (20b) 
150 GeV < m Wu < 225 GeV , (20c) 

where Eq. (20b) is required for only one of the two 6's in the tbj process to be inclusive 3 , 
and all the cuts are in correspondence to the detector level selection criteria stated below. 

Associated tW production is entirely omitted for the time being, because modelling this 
process within its detector acceptance window while at the same time remaining inclusive 
with respect to the other processes is highly nontrivial, and only marginally affects our 
following statements (the main effect being the neglected contamination of the other final 
states at the detector level, which amounts to <10 % in the tb channel and practically 
vanishes in the dominant tj channel) 4 . Clearly, this is still not the fully correct answer at 
detector level, but it should be closer to the truth than entirely neglecting the ^-dependence 
of e, and the consistency of the two approaches can be checked. 

Although it is clear that the object [$ part x a part ] i (g) to be computed is much harder to 
handle than the k's (even more so if the full matrix element response including all off-shell 
and interference effects is to be taken into account), it is basically just a technical issue 
which can be tackled with appropriate Monte Carlo machinery and respective CPU time. 
Therefore, we compare the results of Eq. (18) and Eq. (19), and see if the effort is justified. 

B. Comparison of On-shell and Off-shell Results 

For the measurement of the total cross section of a given final state at the detector level, 
the experimental sensitivity is given in terms of a measure for Aa/a, where estimations 

3 Clearly, for full inclusiveness one would have to entirely drop the distinction of light and b flavors at 
partonic level, but the actual chance of mistagging the light forward jet is negligible once the full event 
topology (cf. tj selection below) is taken into account. 

4 Moreover, the suppression of the huge irreducible tt background in the radiative correction diagram tWb 
still is a topic of vivid discussion (cf. e. g. [18]), a problem which again does not affect the main statements 
of this paper. Still, it is clear that in the end also this channel should be addressed and included in a 
complete study. 
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for total uncertainties are adopted from [36] for consistency, amounting to 20.8 % (tb sel.) 
and 13.5% (tj sel.) for 10fb _1 of LHC data at \fs = 14TeV. Therefore, we separate 
the overall normalization, which is basically given by higher-order SM results for the total 
production cross sections, from the modelling of the normalized LO matrix element response 
as a function of the anomalous coupling set g, i. e. the k coefficients in the on-shell approach 
or, more generally, a function Aa/a(g) = n(g) for each partonic input process i, where 



cf. Eqs. (18) and (19). Since the aim is to accommodate all coupling dependence therein, 
it is fruitful to first set e = 1 and compare this function for the different approaches at the 
partonic level. 

Partonic level 

To be self-consistent, we essentially redo the analysis procedure presented in [36] for the 
on-shell approach, employing Whizard to compute the coefficients of K on as well as produce 
parton level samples, which are then processed with Pythia and Delphes to retrieve e. 
The quadratic coefficients (i. e. ~ gf) are obtained in each production channel, cf. Fig. 1, 
by separately setting each g^ — 1, integrating the total cross section for on-shell single 
top production and finally normalizing to the SM point (V& = 1, Vr = g^^R = 0). The 
interference terms are computed similarly, setting always two couplings to 1 and subtracting 
the quadratic parts from the result before normalizing to the SM. The implementation of 
the vertices and phase space integration has been checked by switching off the pdf's and 
comparing Whizard to analytical results. Using the pdf's and parameter setup quoted 
in [36], Whizard also reproduces the on-shell k coefficients stated there within numerical 
uncertainties. For all further Whizard results, we set m t = 173.1 GeV, mj = 4.2 GeV, 
m w = 80.42 GeV and choose CTEQ6L1 [52] for pdf's. 

In the off-shell approach, the matrix element response is modelled according to Eq. (21b) 
by applying the acceptance cuts $ part , Eq. (20), to the full partonic off-shell matrix elements. 
In this approach, taking all finite width and interference effects into account, it is a priori 




(21a) 



off-shell: < ff (g) 



[$P art x o" part ]. (g) 



(21b) 



[$part x a part] . | gM ' 
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not clear that the function K g (g) obeys a simple polynomial expansion in g, so rather 
than assuming a specific functional form, we use the Whizard machinery to scan the 
entire 4-dimensional parameter space g = (Vl, Vr, g^, gu) (effects of V£ s will be addressed 
separately in Sec. Ill C) within the relevant numerical ranges < Vl < 1.2, —1.2 < Vr < 1.2 
and —0.6 < g^,R < 0.6 (cf. [36]), also including the full dependence of the top width on the 
couplings T t (g). Since T t has already been measured, the most recent bound from D0 being 
T t = 2. 00^43 GeV [53], it is included in our analysis as an additional observable (whose 
sensitivity on g is effectively always superseded by the cross sections). The numerical results 
can then be used to test the validity of the polynomial parametrisation, Eq. (21a), in the 
following way: The normalized matrix element response k may always be expanded as 

K (g) = Y1 Kl (9i) + Yl K2 (&> 9j) + ^^3(9t,9j,9k) + (22) 

i i,j i,j,k 

where the K{ are polynomials in their respective arguments. Offsets kq from squared irre- 
ducible background diagrams could be considered, but are obviously independent of g and 
merely add to the background normalization, so we just subtract them from the scans, while 
keeping all interference effects (appearing as terms linear in the in Ki) for completeness. 
Usually this series terminates after K2, which becomes obvious when applying the narrow 
width approximation, where additional coupling effects cancel in the interplay of the decay 
vertex insertion, width dependence and phase space integration. This basically leads to the 
quadratic form in the on-shell approach. However, in the special case of single tops, produc- 
tion and decay are interrelated via the same set of CC couplings, thus affecting production 
as well as decay distributions, which in combination with the detector acceptance cuts might 
lead to deviations from the on-shell result in some regions of the parameter space 5 . 

In order to estimate the size of the moments K3, in our cross section scans we consider 
2-dimensional subplanes (gi,gj) among the anomalous couplings (Vr, gi, gu) for different 
fixed values of Vl- After subtracting all the lowest moments, 

A«(<7) = K(g) - ^2 K i(9k) - K 2 (g k , V L ) , (23) 

k=i,j,VL k=i,j 

within the on-shell picture the remaining contribution An(g) ~ K 2 (gi,gj) should then be 
independent of the value of Vl, otherwise it would by definition contain some K 3 (gi, gj,Vi)- 

5 Note that this is a qualitative difference to QCD (LHC) or NC (ILC) tt production, where production 
and decay are affected by different sets of anomalous couplings. 
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FIG. 2: A/c|y L= i.2 — A K |vi=o.2 in the ij process (left) and ibj process (right), based on full matrix 
elements (top) and resonant diagrams (bottom). 

For illustration, we choose the coupling subspace (gi,gj) = (Vr, gi), because it is one of the 
dominant interference contributions to all production processes, and evaluate A/«|y i= i.2 — 
Ak|v l= o.2- The resulting picture is mixed: while in the s channel the result is practically 
all over the Vr-ql plane, in the t channel process ij, which plays a central role because 
of its comparably large cross section, this difference amounts to ~ —0.2 at Vr ~ ±1 and 
gi ~ ±0.5 along the interference direction Vr ~ 2gi (cf. Fig. 2), which is of the same size 
as the respective on-shell contribution ~ —0.5 x Vr x g L , (the same is true for the V^gR 
interference in the tj channel). For comparison, we repeat the whole procedure selecting 
only the resonant single top diagrams for the scan (including the full top width dependence 
on g), finding that background interference only plays a minor role. 

Furthermore, note that there are other interference directions also showing substantial 
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FIG. 3: k scan results and quadratic fits of resonant matrix elements inside the acceptance region 
along the coupling direction Vr ~ 2g^, for Vl = 0.6, 1.0 and 1.4 (bottom to top) in the ij channel 
(left) and the ibj channel (right). Grey dashed lines indicate the on-shell k function. 

effects, e. g. in the Ql-Qr plane of the s channel, but the VrQl interference is the most 
interesting one because it is large in all channels, and respective bounds are expected to 
remain rather weak also from other experiments along the considered direction Vr ~ 2ql [36]. 
All in all, these numbers indicate strongly that the interference and off-shell effects contained 
in the full matrix elements might become important at the level of the coupling limits to 
be expected experimentally, at least in a stand-alone single top cross section study at the 
LHC, and should be checked in any case. 

In order to further quantify this effect and exclude potential artefacts from unphysical 
regions in the parameter space, we now systematically scan the VrQl interference along the 
direction Vr ~ 2#£ as a function of Vl, including resonant off-shell diagrams for anomalous 
single top production and decay. We observe that the scans do show a quadratic dependence 
on Vr and ql to a very good approximation at any value of Vl (cf. Fig. 3), so the respective 
coefficients Ky R ~ Vr, k 9l ~ g\ and Kv R g L ~ Vr9l can be extracted from quadratic fits 
along the axes Vr, Ql and Vr = 2ql, for each fixed value of Vl and gR = 0. This is done for 
the full phase space integration as well as for the acceptance region denned in Eq. (20). As 
illustrated in Fig. 4 for the t channel processes, when integrating over the full phase space 
the deviations from the on-shell result are very small as expected from the narrow width 
approximation, whereas application of the acceptance cuts leads to substantially different 
fit coefficients which also depend on the value of Vl- 

Particularly this dependence of the magnitude of the VrQl interference on Vl is an ex- 
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FIG. 4: Numerical results for various k's in the processes tj (left) and ibj (right), as extracted 
from quadratic fits to the normalized cross sections. The dashed line indicates the value of the 
on-shell k in each case. 

ample for the breakdown of the quadratic on-shell approach. However, this dependence is 
rather small in Ky R and k 9l which come with squares of the respective couplings and there- 
fore dominate the sensitivity of a given process to these couplings, so it might still be a 
good approach to expand the full k, q as a quadratic form in small anomalous deviations 
from the SM point Vl = 1, namely by choosing Vl = 1 as an origin. As explained before, 
the numerical values are extracted from quadratic fits to 1-dimensional coupling scans (in- 
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FIG. 5: The differences \K Q g — K on \ (left) and \n g — Kg t | (right) in various coupling planes of 
the tj process. Note that the heat scales on the left are of the size of the expected experimental 
sensitivity to the t channel (~ 13%), and significantly decrease on the right. 
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FIG. 6: The differences \k b — K> on \ (left) and \k q q — Kg t | (right) in various coupling planes of the 
ibj process (note the heat scales, cf. discussion in Fig. 5). 
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eluding acceptance cuts) in all directions of g, and all interference directions g%gj. A major 
difference to the original on-shell quadratic form is the appearance of large linear terms in 
the couplings, which now encode the interference with the SM. In Figs. 5 and 6, we compare 
the different quadratic parametrisations, namely K on from the on-shell approach and K&t 
inferred from the fits, to the full off-shell response k s, illustrating two statements: Firstly, 
n on significantly deviates from k q q in various parts of the parameter space relevant to the 
analysis, particularly in the Ql-Qr and Vr-Ql planes. Secondly, although still neglecting the 
higher coupling dependences, K fit inferred from quadratic fits to the full scan does show a 
significantly improved agreement with the full scan while still being fast and efficient. This is 
further illustrated in Fig. 7 and 8 showing ±lcr sensitivity contours around the SM strength 
k — \ for various anomalous coupling combinations and production processes: Especially 
when the momentum-dependent couplings gi t R are involved, there are remarkable shifts of 
the contours when going from the on-shell to the full off-shell approach, but generally these 
effects are modelled very well by the adapted quadratic parametrisation Kfa(g), while the 
machine cost reduces from a 4D scan over g to a set of ID scans along all axes gi and inter- 
ference directions gigj for each input process. We now go on to the detector level to quantify 
the impact of these different approaches on exclusion bounds on anomalous couplings from 
combined cross section measurements. 

Detector level 

In order to obtain a detector efficiency matrix in the various approaches, samples of 500 k 
events are produced in each partonic production process, once integrating the full off-shell 
matrix elements over the acceptance region, Eq. (20), and once integrating the resonant 
matrix elements over the full phase space, letting the tops decay off-shell into a b quark, 
a charged lepton and a neutrino, analogously to [36]. All the parton level samples are 
processed with Pythia and Delphes to obtain events at detector level. This is done for 
the SM point as well as the coupling configurations 

A: V L = 1, V r = 0.3, g L = 0.15, 

(24) 

B: V L = 1, g R = 0.024, 

taken from [36] to facilitate comparison. Note that in each case a consistent top 
width Y t (g) is calculated beforehand and the result checked to comply with experimen- 
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FIG. 7: Comparison of the la contours for the various matrix element response functions n on , Kfn 
and k s at parton level, for tj and ij production processes in different coupling planes (setting the 
others to their SM values). 

22 



0.6 ^ 
0.4 - 
0.2 - 



tb j channel 
as nn — on -shell 
> ■ — quadr. fit 
— full ME 

-0.2 - 



-0.4 



-0.6 t 




0.6 0.8 1.0 1.2 




0.8 1.0 



tb j channel 
■■■ on -shell 

— quadr. fit 

— full ME 

1.2 



0.6 
0.4 
0.2 
£3 0.0 
-0.2 
-0.4 
-0.6 




t bj channel 
■■■ on -shell 

— quadr. fit 

— full ME 



0.6 



0.8 1.0 

V L 



1.2 



0.6 F 
0.4 
0.2 
£ 0.0 
-0.2 
-0.4 
-0.6 









t bj channel 




— on -shell 




— quadr. fit 




— full ME 







-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

V R 



0.6 








0.6 




0.4 









0.4 




0.2 








0.2 




0.0 








5 o.o 




-0.2 








-0.2 




-0.4 
-0.6 






tb j channel 

— on -shell - 

— quadr. fit 

— full ME 


-0.4 
-0.6 


t bj channel 
■■■ on -shell 

— quadr. fit 

— full ME 



-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

gh 



-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

gL 



FIG. 8: Comparison of the la contours for the various matrix element response functions K on , Kfi t 
and K fi at parton level, for tbj and ibj processes in different coupling planes (setting the others to 
their SM values). 
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tal constraints [53]). Again for comparison reasons, we also adopt the final state selection 
cuts stated in [36] which define the various components of e: apart from requiring an isolated 
lepton (that is, e or /i) with p T > 25 GeV and missing transverse energy $ T > 25 GeV, the 
selection criteria for the three final state signatures are, respectively 

1. for tb selection: exactly two b tagged jets (assuming a tagging efficiency of 0.6) with 
Pt > 30 GeV, and neither central nor forward light jets with p? > 15 GeV. In ad- 
dition, the top momentum p t is reconstructed from one of the 6's together with the 
charged lepton and $ T (to be identified with the neutrino pr), by applying the on- 
shell constraint (pe + p u ) 2 = m^y and picking the smaller of the two solutions for the 
longitudinal component of p v . Finally, the resulting top mass must lie between 150 
and 225 GeV. 

2. for tj selection: at least one b jet with pt > 30 GeV (one of them reconstructing p t 
together with the leptons as explained above), one light forward jet with px > 50 GeV 
and 2.5 < \rj\ < 5 and no more than one additional light central jet, which may have 
p T < 30 GeV only. 

By applying every final state selection to each of the 500 k event samples corresponding to 
the partonic input processes and averaging over lepton flavors and charge states, we find 
for the samples from resonant diagrams integrated over the full phase space an efficiency 
matrix e (in %) at the SM point: 





tb tj tbj 




tbsel. 


0.658(6) 0.040(1) 0.051(2) 


(25) 


tj sel. 


0.165(3) 0.647(6) 0.531(5). 





For the off-shell approach, we run the selection criteria described above on the detector level 
samples from full matrix elements integrated over the acceptance region, inferring e' (in %) 





tb 


tj tbj 


tb sel. 


1.28(1) 


0.039(1) 0.031(1) 


tj sel. 


0.282(4) 


1.52(1) 1.023(7) 



at the SM point. 

Before moving on to detector level coupling limits, the detector efficiency matrices can be 
compared between the various coupling points, Eq. (24), to test the assumption of a constant 
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detector response. While the small admixture of gn in point B only has a negligible effect 
on the efficiencies, we observe that there are sizeable changes in the efficiencies when going 
from the SM point to point A, amounting to some 15 % (8 %) of the original values for tb 
(tj) selection in the samples with full phase space integration at parton level. However, 
when going close to the detector acceptance region already with the partonic input, this 
dependence is reduced to ~ 6% (4%), thus improving on another source of systematic 
uncertainty. 

Taking approximate NNLO on-shell s and t channel production cross sections from the 
literature [54, 55] (multiplied by a partonic acceptance efficiency corresponding to the cuts 
in Eq. (20) in the off-shell approach) to normalize the SM reference point for each input 
process, we now have all ingredients at hand to derive limits on g from a set of cross 
section measurements, modelled by Eqs. (18) and (19), and compare the results. In the 
t channel, the matrix element response for the detector level analysis is modelled using 
only tbj and ibj processes for simplicity, and because it was argued that the corresponding 
distributions already describe the proper NLO behaviour rather well [56, 57]. As already 
anticipated from the la contours in Fig. 7 and 8, while effects on the bounds on Vl,r remain 
small in general, the largest differences are found when the momentum dependent couplings 
gL,R are involved, particularly in the girQR plane illustrated in Fig. 9. In this case, when 
single channels and charge states are considered separately, the different approaches tend to 
produce very different exclusion bounds. Fig. 9 might also suggest that after combining all 
channels and considering the ratio R(i/t) of cross sections for t and t production in the t 
channel as an additional observable (tentatively assuming 2 % statistical and 3 % systematic 
uncertainty as estimated for 10 fb _1 at \fs = 14 TeV in [36] , in the absence of a more detailed 
experimental assessment), one might end up with the same exclusion limits again, but indeed 
this depends heavily on the total uncertainty of R in the actual experiment (cf. "-R2" in 
Fig. 9). In any case, when the aim is to properly understand and quantify the sensitivities to 
anomalous couplings of the various final states separately, going from the on-shell approach 
to full matrix element responses inside the selection acceptance region produces considerable 
effects that should not be neglected. In that respect, the adapted quadratic parametrisation 
introduced above, employing quadratic fits to off-shell scans inside the acceptance region, 
represents a very good approximation to the full off-shell parameter scan (cf. Fig. 5-8). 
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FIG. 9: Combined la limits on gi and gR (Vl = 1, Vr = 0) from single final states (top), 
combined final states (center left), and including the observable R(i/t) (center right and bottom). 
"i?2" denotes a factor 2 on the experimental resolution of R. 
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C. Pinning Down the Off-shell Coupling 

After discussing the technical issue of modelling the LO matrix element response to 
anomalous top couplings at an experimentally relevant acceptance level, and validating 
an adapted quadratic parametrisation which simultaneously meets the demands of machine 
efficiency and good agreement with the full off-shell coupling scan in the previous section, the 
closing section of the article is devoted to the application of the new approach to a physical 
issue, namely a possible admixture of the additional anomalous coupling V£ introduced in 
Sec. II B, to the single top cross sections. 

Considering the experimental sensitivities to the anomalous couplings of the total cross 
sections stated above for the LHC (which are already dominated by systematics) , it is 
clear that a stand-alone study of single top cross sections alone will never provide the most 
stringent bounds on the complete parameter space of anomalous CC couplings, including 
VI or not. Therefore, rather than just adding another direction to g, the focus shall be 
directed here to those regions of the parameter space where single top cross sections actually 
become the crucial inputs to the combined limits. 

More explicitly, the top decay observables (mostly related to the charged lepton distri- 
bution) are very sensitive to anomalous W helicity fractions generated by Vr, and gu 
(cf. [36, 39, 40, 58]). For example, the limit \gn\ ^ 0.024 stated in [36] for our LHC reference 
point (10fb _1 at 14TeV) is more than an order of magnitude below the sensitivity of the 
cross sections, so we may as well set g^ = for our purposes. On the other hand, the large 
interference among Vr and gi leads to rather poor bounds \Vr\ ^ 0.3 resp. \gi\ ^5 0.15 as 
long as they are fine-tuned to Vr ~ 2gL- Finally, since decay observables basically measure 
helicity fractions, they are neither sensitive to the overall vertex normalization nor to the 
admixture of V£ to the left-handed vector part. This is where the cross sections come 
into play, delivering the most stringent direct constraints. In Fig. 10, we therefore present 
combined limits on Vj, and V£ s from single top cross sections, both setting Vr — gz, = 
as well as varying over —0.3 < Vr = 2gi < 0.3. The very different sensitivities of the two 
final states greatly help in the combined limit: the s channel is very sensitive along V£ G due 
to the kinematics, whereas the t channel basically cuts the substantial interference in the s 
channel along Vl- Still, the resulting limit on Vl deteriorates from 0.9 < Vl < 1.1 {V£ s = 0) 
to 0.82 < Vl < 1-1 iVl s varied). Naturally, projecting over the remaining freedom in Vr 
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FIG. 10: Combined la contours in the Vl~V£ s plane, setting Vr = gi = (top) or projecting over 
the direction Vr = 2gi (bottom). 



and gi instead of switching them off further relaxes the combined limits to 0.68 <Vl< 1.1. 
Fig. 11 displays combined bounds in the Vr-Ql plane, switching off resp. varying over V£ ff . 

In the long run, it is perfectly clear that this ambiguity among Vj, and V£ s remaining 
in the single top cross sections can be further resolved, namely by examining differential 
distributions, since Vf scales very differently with the partonic y/s than Vl- (In fact, V£ s 
behaves like, or parametrises, a heavy off-resonant new degree of freedom, cf. e. g. [59].) 
However, this issue will have to be tackled in the s channel where the momentum of the 
W propagator producing the top becomes timelike. Sensitive observables would obviously 
be the total invariant mass m t b of the final state or the pseudorapidity rjb of the hard b jet 
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produced along with the single top. However, such a study is experimentally challenging, 
since it requires a very good isolation of the tiny s channel signal from the huge t channel 
contamination, whereas at present this signal has not even been established yet individually 
at the LHC (the most recent search being [60]). Hence, one should stay careful when stating 
limits on Vl from measurements of the overall size of V t b until its kinematic behaviour is 
further clarified experimentally. 
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IV. CONCLUSIONS 



In this article, we have revisited the model-independent parametrisation of anomalous 
top couplings to SM gauge bosons within the effective operator approach, paying special 
attention to the charged-current sector and its phenomenological implications at current 
hadron colliders. More explicitly, addressing the minimal fully general set of anomalous 
trilinear tbW couplings coming from dimension six effective operators, there is a controversy 
regarding the meaning of "fully general" , namely whether an off-shell interaction contained 
in the original operator basis should be dropped because it turns out to be related to four 
fermion contact interactions after application of the equations of motion, or not. While 
dropping it and sticking to the usual coupling basis (Vl,Vr, <?l, Qr) naturally simplifies the 
analysis, there are good arguments to include it: Since it emerges from the minimal gauge 
invariant operator basis that also generates the trilinear couplings, the respective coupling 
size V£ s is related to the other ones by the underlying operator coefficients. (For example, 
an experimental limit on 5Vl is ambiguous in the context of effective operator coefficients, 
requiring knowledge about either 8V£ S or the NC sector to be resolved.) For the same 
reason, the coupling basis (Vl, Vr, gi, 9r, V£ s ) parametrises the complete set of anomalous 
diagrams which interfere with the SM diagram in a minimal way, so including it is also 
consistent at the level of matrix elements. 

In the phenomenological part, the dependence of single top cross sections on anomalous 
tbW couplings in s and t channel production is examined, stressing the fact that the couplings 
do not only affect the total cross sections but also final state distributions, which determine 
the selection efficiencies within the detector acceptance region. While these effects are 
usually considered small, working with constant detector efficiencies and modelling the whole 
coupling dependence on the basis of on-shell production amplitudes, we use the Whizard 
machinery to scan the full off-shell matrix element dependence on the couplings inside the 
acceptance window defined by the final state selection cuts. Comparing to the on-shell 
approach, one finds considerable deviations in some regions of the parameter space, especially 
where the momentum-dependent couplings ql,r are involved, affecting the sensitivities of 
the various production channels to those couplings and therefore also the limits derived 
from the experiment. Finally, an adapted polynomial approach of the coupling dependence 
is discussed, which is based on quadratic fits to the full off-shell matrix element response 
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including detector acceptance, and turns out to parametrise the full scan result rather well 
while still being fast and efficient. However, it is also stated as a result of the present study 
that the theoretical modelling of the coupling dependence should be adapted as closely as 
possible to a given experimental analysis with defined selection criteria to minimize the 
systematic uncertainty of the derived limits. 

The study concludes with a short discussion of the influence of top decay observables 
on combined coupling limits, and the regions of the parameter space where single top cross 
sections still provide the crucial input to the bounds, namely the overall tbW vertex nor- 
malization and the interference direction Vr ~ 2gi. In this respect, we address the impact 
of including V£ s in the coupling basis, and briefly point out the possibilities to resolve the 
ambiguity between Vl and V£ 5 experimentally, using kinematic distributions in the s chan- 
nel. 
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